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Abstract: We consider a finite horizon linear discrete time varying system whose input is a random noise 
with an imprecisely known probability law. The statistical uncertainty is described by a nonnegative 
parameter a which constrains the anisotropy of the noise as an entropy theoretic measure of deviation 
of the actual noise distribution from Gaussian white noise laws with scalar covariance matrices. The 
worst-case disturbance attenuation capabilities of the system with respect to the statistically uncertain 
random inputs are quantified by the a-anisotropic norm which is a constrained operator norm of the 
system. We establish an anisotropic norm bounded real lemma which provides a state-space criterion 
for the a-anisotropic norm of the system not to exceed a given threshold. The criterion is organized as 
an inequality on the determinants of matrices associated with a difference Riccati equation and extends 
the Bounded Real Lemma of the J£L-control theory. We also provide a necessary background on the 
anisotropy-based robust performance analysis. 
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1. INTRODUCTION 

The statistical uncertainty, present in random disturbances as 
a discrepancy between the imprecisely known true probability 
distribution of the noise and its nominal model, may corrupt 
the expected performance of a stochastic control system if the 
controller design is oriented at a specific probability law of the 
disturbance. Such uncertainties result not only from the lack 
of prior knowledge of the actual noise statistics, but also from 
the inherent variability of the environment where the control 
system operates. 

The robustness in stochastic control can therefore be achieved 
by explicitly incorporating different scenarios of the noise dis- 
tribution into a single performance index to be optimized. The 
degree of robustness depends on the "size" of the uncertainty 
used in the controller design. The statistical uncertainty can be 
measured in entropy theoretic terms and the robust performance 
index can be chosen so as to quantify the worst-case distur- 
bance attenuation capabilities of the system. It is this combi- 
nation of approaches that underlies the anisotropy-based theory 
of stochastic robust control which was initiated about sixteen 
years ago at the interface of the entropy /information and robust 
control theories in a series of papers (Semyonov et al. (1994); 
Vladimirov et al. (1995a,b, 1996a,b, 1999)). This theory em- 
ploys the anisotropy functional as an entropy theoretic measure 
of deviation of the unknown actual noise distribution from the 
family of Gaussian white noise laws with scalar covariance 



EAM and APK are supported by the Russian Foundation for Basic Research 
grant 08-08-00567. IGV is supported by the Australian Research Council. 
' We deeply regret the untimely loss of Dr Eugene Maximov who passed away 
on 26 July 2010. His talents and ideas will never be forgotten. 



matrices. Accordingly, the role of a robust performance index 
is played by the a-anisotropic norm |||F||| a of a system F which 
is defined as the largest ratio of the root mean square (RMS) 
value of the output of the system to that of the input, provided 
that the anisotropy of the input disturbance does not exceed a 
given nonnegative parameter a. Thus, the input anisotropy level 
a is the size of the statistical uncertainty, and the a-anisotropic 
norm of the system |||F||| a is the worst-case RMS gain which, in 
the framework of the disturbance attenuation paradigm, is to be 
minimized. 

An important property of the a-anisotropic norm is that it 
coincides with a rescaled Frobenius (or J^) norm of the system 
for a = and converges to the induced (or ML) norm as a — > 
+°°. Therefore, ||| • ||| is an anisotropy-constrained stochastic 
version of the induced norm of the system which occupies a 
unifying intermediate position between the J#2 and J^o-norms 
utilised as performance criteria in the linear quadratic Gaussian 
(LQG) (Kwakernaak & Sivan (1972)) and J£L-control theories 
(Doyle et al. (1989)). 

In its original infinite horizon time invariant setting, the 
anisotropy-based theory employed the anisotropy production 
rate per time step in a stationary Gaussian random sequence. 
The mean anisotropy has useful links with the condition num- 
ber of the covariance matrix and the transient time in the se- 
quence, thus describing the amount of spatial non-roundness 
and temporal correlation in it. These connections have recently 
been revisited in (Kurdyukov & Vladimirov (2008)). 

At the performance analysis level, the anisotropy-based theory 
was developed for time invariant systems in (Vladimirov et al. 
(1996a)), where equations were obtained for computing the a- 
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anisotropic norm in state space. An extended exposition of that 
work can be found in (Diamond et al. (2001)) and a general- 
ization to finite horizon time varying systems is provided by 
(Vladimirov et al. (2006)). 

A state-space solution to the anisotropy-based optimal control 
problem (which seeks an internally stabilizing controller to 
minimize the a-anisotropic norm of the closed-loop system) 
was obtained in (Vladimirov et al. (1996b)). The solution was 
found as a saddle point in the stochastic minimax problem. 
The anisotropy-based theory therefore offers tools both for the 
quantitative description of uncertainty and robust performance 
analysis on the one hand, and control design on the other. 

The present paper is concerned with the anisotropy-based 
robust performance analysis of linear discrete time varying 
(LDTV) systems in state space. The procedure developed in 
(Vladimirov et al. (2006)) for computing the a-anisotropic 
norm of such a system over a finite time horizon involves 
the solution of three coupled equations: a backward difference 
Riccati equation, a forward difference Lyapunov equation, and 
an algebraic equation for the determinants of matrices associ- 
ated with the previous equations. Its practical implementation 
is complicated by the opposite time ordering of the coupled 
difference equations. 

Here, we develop an alternative state-space criterion for the 
a-anisotropic norm to be bounded by a given threshold. This 
allows the above issue to be overcome by eliminating the 
Lyapunov equation, replacing the backward Riccati equation 
by a forward Riccati equation, and replacing the algebraic 
equation by an appropriately modified inequality. The resulting 
Anisotropic Norm Bounded Real Lemma (ANBRL) is organ- 
ised as an inequality on matrices associated with a forward 
difference Riccati equation. In addition to a significant simplifi- 
cation of the previously developed anisotropy-based robust per- 
formance analysis (which can now be carried out recursively in 
time), ANBRL also provides an extension of the Bounded Real 
Lemma from the JtfL -control theory to the uncertain stochastic 
setting with finite horizon time varying dynamics. An infinite- 
horizon version of ANBRL for time invariant systems, which 
involve algebraic equations free from the opposite time order- 
ing issue, is presented in (Kurdyukov et al. (2010)). 

Another approach to robust control in stochastic systems, using 
the relative entropy to describe statistical uncertainty, can be 
found in (Petersen et al. (2000); Petersen (2006); Ugrinovskii & 
Petersen (2002)), where an important role is played by a link be- 
tween a relative entropy duality relation and robust properties of 
risk-sensitive controllers (Dupuis et al. (2000)) which minimize 
the expected-exponential-of-quadratic functional. Although the 
ideas of entropy-constrained induced norms and associated 
stochastic minimax find further development in the control lit- 
erature (Charalambous & Rezaei (2007)), the anisotropy-based 
theory of stochastic robust control remains largely unnoticed. It 
is partly for this reason that the main result of the present paper 
is preceded by the background material to assist the readers to 
build awareness of the anisotropy-based approach. 

The paper is organised as follows. Section 2 specifies the class 
of systems being considered. Sections 3 and 4 provide the 
necessary background on the anisotropy of random vectors and 
the a-anisotropic norm of matrices. Section 5 establishes the 
Anisotropic Norm Bounded Real Lemma. Its connection with 
the Bounded Real Lemma in the limit a -> +°° is discussed 



in Section 6. Section 7 gives concluding remarks. Appendix 
provides a subsidiary state space criterion of outerness. 

2. CLASS OF SYSTEMS BEING CONSIDERED 

We consider a linear discrete time varying (LDTV) system F 
on a bounded time interval [0,N] . Its n-dimensional state x k and 
r-dimensional output z k at time k are governed by the equations 

x k +\=A k x k +B k w k , (1) 
z k = C k x k +D k w k , (2) 

with initial condition xq = 0, which are driven by an m- 
dimensional input w k . Here, A k , B k , C k , D k are appropriately 
dimensioned real matrices which are assumed to be known 
functions of time k. The state-space equations (l)-(2) are writ- 
ten as 



(3) 



where we have also shown the dimensions. For any two mo- 
ments of time s ^ f, the values of the input and output signals 
W and Z on the interval [s,t] are assembled into the column- 
vectors 

W s ,~[wl...wJ} J , Z„ f :=[zJ,...,z f T ] T . 

Since the state of the system is zero-initialized, then Zo :r = 
Fq-jWqj, where F s:t is a block lower triangular matrix with 
(rxra)-blocks fj k given by 

(Cj<t> hk+l B k for j>k 
F S :t:=block s<jtk<t (fj k ), f jk =< D k for j = k. (4) 
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otherwise 



(5) 



is the state transition matrix from x k to Xj for j ^ k, with 
^kk = h% the identity matrix of order n. Since the matrix Fq : n 
completely specifies the system F on the time interval [0,N] as a 
linear input-output operator from Wq.n to Zq-n, all the norms of 
F are those of Fq-n- In particular, the finite-horizon counterparts 
of the J?2 and ^C-norms are described by the Frobenius and 
operator norms of Fq :N as 



11*112:= 



■■Omax{Fo:N), (6) 



where <7 max (-) is the largest singular value of a matrix. We 
will be concerned with the a-anisotropic norm of the system F 
which is also understood in terms of the matrix Fo:iV- This norm 
is obtained by modifying the concept of induced norm with the 
aid of an additional constraint on the input which involves the 
entropy theoretic construct of anisotropy. 

3. ANISOTROPY OF RANDOM VECTORS 

The relative entropy (or Kullback-Leibler informational diver- 
gence) (Cover & Thomas (2006)) of a probability measure P 
with respect to another probability measure M on the same 
measurable space is defined by 

D(P||M):=Eln^. 

Here, P is assumed to be absolutely continuous with respect 
to M with density (Radon-Nikodym derivative) dP/dM, and E 
denotes the expectation in the sense of P. The relative entropy 
D(P||M), which is always nonnegative, vanishes only if P = M. 
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In what follows, D(P\\M) will also be written as D(^||tj) or 
D(/||g) if the probability measures P and M are distributions 
of random vectors E, and 77 or are specified by their probability 
density functions (PDFs) / and g with respect to a common 
measure. 

For any X > 0, we denote by p e x the i?-variate Gaussian PDF 
with zero mean and scalar covariance matrix Xlf 

Pu {w) = {2nX)- l l 2 z-\ w \ 2 l {2X \ weR e . (7) 

Let W be a square integrable absolutely continuous random 
vector with values in W and PDF /. Its relative entropy with 
respect to the Gaussian probability law (7) is computed as 



D(/||p a )=Eln 



f(W) 
Ptx(W) 



= £ln(2*A) + ffl^-h(n (8) 



where 



h(W) := -Eln/(W) = - / /(w)ln/(w)dv 



is the differential entropy (Cover & Thomas (2006)) of W. 
For what follows, the class of square integrable absolutely 
continuous Revalued random vectors is denoted by \J 2 - 

Definition 1. (Vladimirov et al. (1995a, 2006)) The anisotropy 
A(W) of a random vector W e \J 2 is defined as the minimum 
relative entropy (8) of its PDF / with respect to the Gaussian 
PDFs (7) with zero mean and scalar covariance matrices: 

A(W) := inf D(/|| m ) = | ln ^ggg) _ h(w) , (9) 

A similar construct to the rightmost expression in (9) was 
considered for scalar random variables in a context of time 
series prediction in (Bernhard, 1998, Definition 4 on p. 291 1). 

The minimum with respect to X in (9) is achieved at X = 
E(\W\ 2 )/£. The corresponding "nearest" zero-mean Gaus- 
sian random vector W* has the covariance matrix cov(Wh.) = 
E(|W| 2 )/(j/^, and its differential entropy coincides with the 
first term on the right-hand side of (9), that is, h(W») = 
nn(27teE(\W\ 2 )/e)/2. The class of Revalued Gaussian ran- 
dom vectors W with zero mean and a given nonsingular covari- 
ance matrix E will be written as G e (L) (it is a subclass of L|). 
Their PDF is 



p(w) = (2a;)- £ / 2 (detE 



"l/2 f 



where ||v||m := ^/Tr(v T Mv) is the Euclidean (semi-) norm of a 
vector v weighted by a positive (semi-) definite matrix M. 

Lemma 1. (Vladimirov et al. (1995a, 2006)) 

(a) The anisotropy A(W), defined by (9), is invariant under 
rotation and scaling of W, that is, A(XUW) = A(W) for 
any orthogonal matrix U € M? x * and any Aet \ {0}; 

(b) The anisotropy of a random vector W € h f 2 with a given 

matrix of second moments E(WW T ) = E satisfies 

1 £E 
A(W) > --lndet — . 
v ; 2 TrE 

This inequality holds as an equality if and only if W is 

Gaussian with zero mean and covariance matrix cov(W) = 

E; 

(c) For any random vector W € Lj, its anisotropy A(W) is 
always nonnegative and vanishes only if W is Gaussian 



distributed with zero mean and scalar covariance matrix 
(that is, cov(W) = Xlt for some X > 0). 

Lemma 1(a) shows that A(W) quantifies the rotational non- 
invariance of the PDF of W. This property originally motivated 
the term "anisotropy" for the functional. In application to 
Gaussian random vectors W, the assertions (b) and (c) of 
the lemma allow A(W) to be interpreted as a measure of 
heteroscedasticity and cross-correlation of the entries of W. 

Furthermore, Lemma 1(b) implies that if an arbitrary random 
vector W <G Lj, with second-moment matrix E := E(WW T ), is 
replaced by a Gaussian vector T with zero mean E(r) =0 and 
covariance matrix cov(r) = E, then the transition W h-> T is an 
anisotropy-decreasing operation which preserves the second- 
moment matrix, that is, A(r) A(W) and E(rr T ) = E. 

Another important property of the anisotropy functional is its 
superadditivity 

A(W) > A(Wi) + A(W 2 ), W := [Wi T ,W 2 T ] T 
with respect to partitioning the random vector W into subvec- 
tors W\ and W2 (Vladimirov et al., 2006, Lemma 3 on p. 1269). 
This superadditivity is closely related to the asymptotically 
linear growth of the anisotropy for long segments of a stationary 
random sequence that allows the mean anisotropy to be defined 
as the anisotropy production rate per time step (Vladimirov 
etal. (1995a)). 

4. a-ANISOTROPIC NORM OF MATRICES 

Let F € W xi be an arbitrary matrix. We will interpret it 
as a deterministic linear operator whose input is a square 
integrable Revalued random vector W which is considered to 
be a disturbance. While the disturbance attenuation paradigm 
seeks to minimize the magnitude of the output Z := FW, the 
probability distribution of W can be regarded as the strategy of 
a hypothetical player aiming to maximize the root-mean-square 
(RMS) gain of F with respect to W: 



R(F,W) 



E(|Z|- 



(10) 



V E(|W| 2 )' 

Here, the squared Euclidean norm | • | 2 of a vector is interpreted 
as its "energy", so that E(|W| 2 ) and E(|Z| 2 ) describe the 
average energy (or power) of the input and output of the 
operator F, respectively. The denominator E(|W| ) in (10) 
vanishes only in the trivial case, where W = with probability 
one, which is excluded from consideration. 

The map F ^ R(F,W) is a semi-norm in W xf: . It is a norm 
if and only if the matrix of second moments E := E(Wff T ) of 
the random vector W is nonsingular. Indeed, since E(|W| 2 ) = 
TrE and E(|Z| 2 ) = Tr(FT.F T ), the semi-norm properties follow 
from the representation of the RMS gain (10) in terms of 
the Frobenius norm || • H2 as R(F,W) = \\F£l\\2, where Q. := 
y'E/TrE. This also shows that R(F,W) = implies F = if 
and only if Q. is positive definite which is equivalent to the 
positive definiteness of E. 

The RMS gain R(F,W) depends on the matrix F only through 
F T F and never exceeds the induced operator norm ||F||oo. If 
there are no restrictions on the probability distribution of W 
otherthanthe square integrability E(|W| 2 ) < +°°, thenR(F,W) 
can be made arbitrarily close to its upper bound ||F||oo. This 
is achieved by concentrating the distribution of W along the 
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eigen-space of the matrix F T F in M. e associated with its largest 
eigenvalue ||F|| 2 . However, except when the matrix F T F is 
scalar, such distributions are singular with respect to the l- 
dimensional Lebesgue measure and should be considered as 
"non-generic". 

For any random vector W € L|, we quantify "non-genericity" 
of its probability distribution by the anisotropy A(W). Accord- 
ingly, we assume that the disturbance player is constrained 
by the condition A(W) a, where a is a given nonnegative 
parameter. In particular, if a = 0, then by Lemma 1(c), the 
player is allowed to generate only Gaussian random vectors 
W G Ua>o G\XIi) with zero mean and scalar co variance matri- 
ces. With respect to any such W, the RMS gain (10) of the oper- 
ator F becomes a scaled Frobenius norm: R(F,W) = \\F\\ 2 /VI. 

Definition 2. For any a ^ 0, the a-anisotropic norm of a matrix 
F G R sx( is defined as an anisotropy-constrained upper enve- 
lope of the RMS gains (10): 

|||F||| a :=sup{R(F,W): W e h e 2 , A(W) ^ a}. (11) 

This definition closely follows the concept of an induced norm, 
with the only, though essential, difference being the constraint 
A(W) ^ a on the anisotropy (9). It is the latter point where 
the entropy theoretic considerations enter the construct of the 
a-anisotropic norm (11), thus making |||F||| a an anisotropy- 
constrained stochastic version of the induced operator norm 



For any given matrix F € 



JiSXi 



the a-anisotropic norm |||F||| is 



a nondecreasing concave function of a ^ 0, which satisfies 

ll^lb 



\\n\o< 



^ lim F 



(12) 



The rate of convergence of |||F||| to the limiting values 

||F|| 2 /V^ and ||F||oo is investigated in (Vladimirov et al. (1999, 
2006)). The relations (12) show that the a-anisotropic norm 
occupies an intermediate unifying position between the scaled 
Frobenius norm and the induced operator norm. 

Note that b\\F\\ 2 /Vi + (1 - fe)||F||oo, with b := exp(-a), 
also provides an "intermediate norm". However, unlike the 
naive convex combination of the extreme norms, ||| • ||| is an 
anisotropy-constrained operator norm, the very definition (11) 
of which is concerned with the worst-case disturbance attenu- 
ation capabilities of the linear operator (measured by the RMS 
gain (10)) with respect to statistically uncertain random inputs 
(with the uncertainty being measured by the anisotropy (9)) 
and combines both power and entropy concepts in a physically 
sound manner. 

5. ANISOTROPIC NORM BOUNDED REAL LEMMA 

In application to the LDTV system F of Section 2, the a- 
anisotropic norm |||F||| := |||Fo :i v||| a , computed for the matrix 
(4), provides a robust performance index of the system with 
respect to statistically uncertain random disturbances W over 
the time interval [0,N]. In this case, the relations (12), with 
I := m(N + 1 ) the dimension of Wo-n, take the form 

-7== = \n < \\\n\ a < piil = 11*11... (i3) 

^(JV+l) c! ^ +o ° 
where the norms (6) are used. For a time invariant system F, 
its a-anisotropic norm on the interval [0,N] with a := aN tends 
to the a-anisotropic norm of the system as N — > +°°. If the 



anisotropy level a grows sublinearly (a = o(N)) or superlinearly 
(a = 1 /o(l /N)) with the time horizon Af, then |||F||| a converges 
to either of the extreme norms of the time invariant system 
ll-f Ib/v^" or 11*11°°' respectively. 

In (Vladimirov et al. (2006)), computing the anisotropic norm 
of a finite horizon LDTV system in state space was reduced to 
solving three coupled equations: a backward difference Riccati 
equation, an algebraic equation involving the determinants of 
matrices, and a forward difference Lyapunov equation. The 
procedure of the anisotropy-based robust performance analysis 
is complicated by the presence of coupled difference equations 
with opposite time ordering. 

The theorem below provides a state-space criterion for the a- 
anisotropic norm to be bounded by a given threshold y. It turns 
out that the above issue can be overcome by eliminating the 
Lyapunov equation and replacing the algebraic equation by 
an appropriately modified inequality. Moreover, the backward 
Riccati equation can be replaced by a forward Riccati equation. 
By analogy with the Bounded Real Lemma in the ^.-control 
theory, we call the theorem Anisotropic Norm Bounded Real 
Lemma. 

Theorem 2. Let F be an LDTV system with the state-space 
realization (3). Then its a-anisotropic norm on the time interval 
[0,N] satisfies |||F||| a ^ 7 if and only if there exists q > such 
that for the matrices R k €M. nxn , with k = 0,...,N, governed by 
the difference Riccati equation 



Rk+i = A k R k Aj + qB k Bj + M k S k Mj, 



(14) 
(15) 

S k = I r -C k R k C k 1 -qD k D i k , (16) 

with the initial condition Rq = 0, the matrices So, • • • ,Sn are all 
positive definite and satisfy the inequality 



M k = -(A k R k Cj + qB k Dj)S^ 1 



£ IndetS* ^ m(N + l)ln(l - qy 2 ) +2a. 

k=0 



(17) 



Prior to proving the theorem, note that the matrices So, . . . ,Sn, 
defined by (16), are all positive definite if and only if q < 
llFll^ 2 . For any such q, the left-hand side of (17) is nonpositive, 
since S k =4 I m (and so, lndetS,t ^ 0). Hence, any q satisfying the 
specifications of Theorem 2 must also satisfy the inequalities 

-2t' -Tn/m^ —9 a 



y- 2 (l-e- 2a ^ m )^q<y- 2 , a := 



N+l 



(18) 



Here, the ratio a is the anisotropy production rate per time step. 
Therefore, if a significantly exceeds the dimension m of the 
input W, then (18) yields a relatively narrow localization of the 
candidate values for q about y~ 2 . 

Proof. Consider a class W a := {W 0:iV € ^[ ■ a (Wo-.n) «S a} of 
square integrable absolutely continuous random inputs to the 
system F on the time interval [0,N] with the anisotropy (9) 
bounded by a, where I := m(N + 1). By applying (10) and 
(1 1) to the matrix Fq-n in (4), with which we identify the finite 
horizon system F, it follows that the inequality |||F||| ^ y is 
equivalent to the fulfillment of R(F 0:A r, W 0:N ) = y/Tr (Ml) < y 
for all Wq.n € W a . Here, the RMS gain is completely specified 
by the matrices 

E 
TrE' 

with E := E(Wq-.nWq. n ), so that II can be any positive definite 
matrix of order I with unit trace. Application of Lemma 1(b) 



(19) 
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to Wo.N yields the inequality A(Wq-n) ^ — lndet(£Q)/2 which 
becomes an equality if and only if Wq.n is Gaussian distributed 
with zero mean and covariance matrix cov(Wojv) = all for 
some a > 0. Thus, the minimum anisotropy of the disturbance 
Wo-.N, required to achieve a given value 7 of the RMS gain of 
the system, is 



min A(Wo at) : 



1 



max 



2 n^0:Trn=l,Tr(An)^ 



lndet(m) 



(20) 

and is delivered by zero mean Gaussian random vectors Wo-.n 
with covariance matrices proportional to 



y(q):=(I e -qAy 



(21) 



with ^ q < \\F\\„ 2 . Here, assuming that the matrix A is not 
scalar (the trivial case is excluded from consideration), the 
associated functions 



^(1) '■= -r-lndet ^ 
w 2 Try(q) 



jK(q) := 



/ Tr(A^(g)) 
1r,y(q) 



(22) 



(23) 



are strictly increasing in q. This allows the a-anisotropic norm 
of the system to be computed as |||F||| a = jV{.e/~ l (a)), where 
is the functional inverse of si . The solution of the con- 
strained optimization problem (20), described above, is ob- 
tained by using the method of Lagrange multipliers, the Frechet 
derivative <9j;lndet£ = L _1 and strict concavity of IndetE on 
the cone of positive definite matrices £ (Horn & Johnson, 
2007, Theorem 7.6.7 on p. 466); see (Diamond et al. (2001); 
Vladimirov et al. (1996a, 2006)) for details. However, instead 
of computing the a-anisotropic norm |||F||| a , we will proceed di- 
rectly to considering the structure of its sublevel set |||F||| a < 7. 
To this end, from (2 1 ), it follows that A = (I e - y (q) - 1 ) jq and 

Tr(A^( g )) = Tr ^ (g) ^ (24) 

q 

which, in combination with (23), implies that 



Tvy(q) = 



(25) 



Substitution of the last identity into (22) yields 

£f(q)=Ql(q,^(q)), (26) 

where 

2% 7) := \ lndetft - qA) - | ln(l - qf). (27) 

Since — ln(l — qf) is monotonically increasing in 76 [0, 1/ y/q), 
then so is 21(9, 7). In particular, by the strict monotonicity of the 
function Jf , 

z/(q)^Vl(q,Y) for 9 ^V" 1 (7). (28) 
Although the computation of the a-anisotropic norm |||F||| a 
requires both functions srf and ,jV from (22) and (23), the 
function 2l(-,7) contains all the information about the system 
F needed to decide whether it satisfies |||F||| a ^ 7. This is based 
on the property that the function 21(9, 7) achieves its maximum 
at the point q = jY~ x {y) where, in view of (26), it coincides 
with the function stf : 

max 21(9,7) = %(^~\y),Y) = si {JY~ X {y)\ (29) 

0< 9 <||F||^ 2 



as shown in Fig. 1. Before proving this property, we note 
that the inequality ^ 7 is equivalent to si{jV~ l (7)) ^ 




Fig. 1. Typical graphs of si{q) and 21(9,7) from (22) and 

(27) as functions of q. The graphs intersect each other, 
and 21(9,7) attains its maximum, at q = The 
maximum value si (7)) is the minimum anisotropy 
of the input disturbance required to achieve the level 7 for 
the RMS gain of the system in (20). The mutual position 
of the graphs for q ^ jY ~ 1 (7) is qualitatively described by 

(28) . 

a. On the other hand, (29) implies that si(^V~ x {y)) ^ a is 
equivalent to the existence of 96 [0, ||F||~ 2 ) such that 21 (a, 7) ^ 
a. Therefore, 

2%, 7) ^ a for some q •$=>• |||F||| a ^ 7. (30) 

Now, to prove (29), we differentiate the function 21 from (27) in 
its first argument: 



7) =^ 9 lndet(/< - qA) + 



if 



2(1 -qf) 
if 



= - 2 MA-y(q)) + 2{1 _ qf) 



if 



2(l-9^( 9 )2) 2(1 _ qf) 

2{\-qf){\-q^{qfy K ' 

where d q (-) is the partial derivative with respect to q, and use is 
made of (21), (24), (25). From (31) and the strict monotonicity 
of JY , it follows that d q %L(q, 7) is positive for q < ,yV~ x (y) and 
negative for q > ,jV~ x (7), thereby establishing (29). 

We will now express the condition 21(9, Y) ^ a on the function 
(27) in terms of the state-space dynamics of the system (1)- 
(2) which has not been used yet. Recall that, for any con- 
formable matrices M and U, the spectra of MU and UM can 
differ from each other only by zeros (Horn & Johnson, 2007, 
Theorem 1.3.20 on p. 53). Hence, if we change the order in 
which Fq. n and Fq-n are multiplied in the definition of the matrix 
A in (19), then it follows that the spectrum of I( — qA differs 
from that of / r (jv+i) — qFo-JfFQ-N om Y by ones , an d so 

det(I e - 9A) = det(/ r(Af+1) - qF : N F ( j N ). (32) 

Now, the latter matrix can be factorized for any q < ||F||~ 2 as 

h(N+l) ~ iFo:nFo;n = Ho-.nHq-n, (33) 
where the matrix Ho :N represents an LDTV system H, whose 
input and output are both r-dimensional, on the time interval 
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[0,N]. The factorization (33) is equivalent to that an ancillary 
system 

W:=[y/qF IT], (34) 
with an (m + r) -dimensional input and r-dimensional output, 
is outer on the time interval [0,N], that is, ^orA^o/v = I r (N+i)- 
This property means that *P transforms an (m + r)-dimensional 
Gaussian white noise sequence (with zero mean and identity 
covariance matrix) at the input into an r-dimensional sequence 
with the same properties at the output. The system H can be 
found in the form 



H 



A. 




c. 


Vs. 



(35) 



where the matrices Sq,...,Sn are all positive definite, so that, 
in view of (3), the state-space realization of VP in (34) is 



A. 



C 



y^B. M.y/S. 



(36) 



Application of the state-space criterion of outerness from Ap- 
pendix to (36) yields the Riccati equation (14)— (16), where 
R ,...,Rn are the controllability gramians of the system \P. 
Since (35) implies that Hq-m is a block lower triangular matrix , 
that is, 



0:N : 



with the blocks t/Sq, \fS~N over the main diagonal, then (32) 
and (33) yield 

N 

detfa-qA) = (deti/ :iv) 2 = n det5it - 

k=0 

By substituting this representation into (27), it follows that the 
inequality y) a in (30) is equivalent to (17). ■ 

6. INFINITE ANISOTROPY LIMIT 

If the anisotropy level increases unboundedly, a — > +°°, then the 
localization (18), which follows from the inequality (17), yields 
q — )• y~ 2 . In this case, the Riccati equation (14)— (16) takes the 
form 



T 

R k+l = A k R k A k 



-M k S k M T kl 



(37) 

M k = - {A k R k q + y A B k D k )S" 1 , (38) 
S k = I r -C k R k Cj-y- 2 D k D T k , (39) 

well-known in the context of JC-suboptimal controllers. This 
is closely related to the convergence \im a ^ +m \\\F\\\ a = ||F||„ in 
(13), whereby the inequality |||F||| a < /"approaches" ||F||<x, ^ y 
for large values of a. Therefore, in the limit, as a — > +°°, 
Theorem 2 reduces to the Bounded Real Lemma, which estab- 
lishes the equivalence between the inequality ||F || „ < y and the 
positive definiteness of the matrices Sq,...,Sn associated with 
(37)-(39). 

7. CONCLUSION 

We have considered a class of finite-dimensional linear discrete 
time varying systems on a bounded time interval subjected to 
input disturbances with an unknown probability law. 

The statistical uncertainty has been quantified using the concept 
of anisotropy as an entropy theoretic measure of deviation of 



the actual noise distribution from nominal Gaussian white noise 
distributions with scalar covariance matrices. 

The associated robust performance index, describing the worst- 
case disturbance attenuation capabilities of the system, is the 
a-anisotropic norm defined as the largest root mean square gain 
of the system with respect to random noises whose anisotropy 
is bounded by a given nonnegative parameter a. 

We have established a state-space criterion for the a-anisotropic 
norm not exceeding a given threshold value. The Anisotropic 
Norm Bounded Real Lemma (ANBRL) is organised as an 
inequality on the determinants of matrices associated with a 
forward difference Riccati equation. 

Apart from substantially simplifying the previously developed 
procedure of anisotropy-based robust performance analysis, 
which is now amenable to a recursive implementation, ANBRL 
includes the Bounded Real Lemma of the .^C-control theory 
as a limiting case, thus extending it to the statistically uncertain 
stochastic setting with time varying dynamics. 
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Appendix. STATE-SPACE CRITERION OF OUTERNESS 

In this section, we assume that the output dimension of the 
system F in (3) does not exceed the input dimension: r ^ m. 
Such a system is said to be outer on the time interval [0,N] 
if Fq;nFq. n = I r (N+iy This is equivalent to the preservation of 
the property of being a Gaussian white noise sequence (with 
zero mean and identity covariance matrix) for signals passing 
through the system. The state-space criterion of outerness, 
given below for completeness of exposition, is similar to that 
in the time invariant case (Gu et al. (1989)) and utilizes the 
controllability and observability gramians 

£=0 j=k 

where <t>j k is the state transition matrix (5). The gramians 
satisfy the difference Lyapunov equations 

Pj+i = AjPjA] + BjBj, Q k = A T k Q k+l A k + C T k C k 

with initial condition Pq = and terminal condition Qn+i = 0. 

Lemma 3. The LDTV system F with the state-space realization 
(3) is outer on the time interval [0,7V] if and only if for every 
k V. 

C k P k C T k + D k Dj = I r , Q k+l (A k P k C T k + B k D T k ) = 0. (A.2) 

Proof. Let the input W to the system F be an m-dimensional 
Gaussian white noise sequence with zero mean and identity 
covariance matrix. Then, in view of (2), the output Z is also 
a zero mean Gaussian random sequence with 

covfe) = C k P k Cj + D k D T k , (A3) 



since x k and w k are independent, with cov(x k ) = P k the A:th con- 
trollability gramian from (A.l). Consider the cross-covariance 
of Zj and z k for j > k > 0. By using the state transition matrix 
(5), it follows that 

7-1 

Zj=Cj{^j,k+\Xk+\+ t I > /\* |/>\n\! • /)/«'/• 

s=k+\ 

Therefore, since w k+ \,...,Wj are independent of x k and w k , 
then 

COV(zj,Zk) = Cj®jf + iCOV(x k+ i,Zk) 

= C J <P jtk+l (A k P k Cj + B k Dj). 

Hence, by recalling the observability gramian from (A.l), it 
follows that 

£ \\coy {zj,Zk)\\l=\\V^{A k PkCl+B k Dl)\\l (A.4) 
j=k+i 

Now, the system F is outer on the time interval [0,7V] if and 
only if cm{zj,Zk) — Sj k I r for all TV ^ j ^ k ^ 0, with 8j k the 
Kronecker delta. In view of (A. 3) and (A.4), the outerness is 
equivalent to (A.2). ■ 



